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Abstract 

We prove that in dimension one the non-real eigenvalues of the non-Hermitian Ander- 



> 

o 
o 

^ ■ son (NHA) model with a selfaveraging potential are regularly spaced. The class of selfav- 

CN ! eraging potentials which we introduce in this paper is very wide and in particular includes 

stationary potentials (with probability one) as well as all quasi-periodic potentials. It should 

• be emphasized that our approach here is much simpler than the one we used before. It al- 
lows us a) to investigate the above mentioned spacings, b) to establish certain properties 

I of the integrated density of states of the Hermitian Anderson models with selfaveraging 

• potentials, and c) to obtain (as a by-product) much simpler proofs of our previous results 
S . concerned with non-real eigenvalues of the NHA model. 

> 

• 1—1 . 

^ : 1 Introduction 



The non-Hermitian Anderson model (NHA model) was introduced by N. Hatano and D. Nelson 
in 1996. It arises in the physics of vortex matter, and in many other contexts, see e.g. 



[Q, [IT], |T2p. This model is described by the following operator 

{H^Lp)k = -e^^k+i + QkVk - e'^ifk-i, l<k<n (1) 
with periodic boundary conditions 

= ifn, (Pl = (fn+l- (2) 

Here g is a real parameter, g > 0. The Hilbert space is ^2(15 with the standard inner product: 
if (f = {(fj}"^-^ and ijj = {ipj}^^^ are two vectors from /2(1, n), then (v?, ip) = J2]=i ^j'^j- 

Hatano and Nelson considered the case when the values qj of the potential are taken as a re- 
alization of a sequence of independent identically distributed random variables. By conducting 
a numerical experiment they discovered a number of remarkable properties both of the spec- 
trum and the eigenf unctions of the operator H^. It turns out that the asymptotic behavior of the 
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eigenvalues depends strongly on the value of the parameter g. To demonstrate this statement 
we present in Fig.[I] results of a similar numerical experiment. These pictures are not so well 
predictable in the following sense. It is a consequence of the Weyl criterion that the spectrum 
of the limiting random operator {n = oo) contains with probability 1 the union of spectra of 
operators with constant potentials, qj = q, for any real q belonging to the support of the random 
variable qi . A very simple calculation shows then that the spectrum of would typically con- 
tain (with probability 1) a two-dimensional subset of the complex plain. (A much more detailed 
description of the spectrum of the limiting operator can be found in However, numerical 
experiments reproduce pictures like that in Fig.|T] with remarkable stability also for large values 
of n (in ^ n = 1000). They clearly show that the eigenvalues of have no tendency to 
spread over any two-dimensional region but rather tend to belong to smooth curves. 
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Figure 1: The eigenvalues of H^, are presented by dots plotted in the complex plane. Here 
n = 50 and {qk} is a fixed realization of independent samples from the uniform distribution on 
[-4, -3] U [3, 4]; the values of g are: (a) g = 0.2, (b) g = 1.1, (c) g = 1.4. 

Our attempt to understand whether this phenomenon would indeed persist as n — oo or 
is due to the fact that n is not large enough stimulated the appearance of two papers [0, ^ 
where the analysis of the spectra of operators for finite but large values of n was carried out. 
We shall now briefly describe some of our results restricting ourselves to the case of bounded 
potentials. Namely it turns out that there are critical values g^^ and gcr depending only on the 
distribution of the random potential and such that: 

1. If < (7 < g then the eigenvalues of are "asymptotically" real (see Theorem ^]2| 
for the exact meaning of this statement). Moreover their limiting distribution does not depend 
on g and is the same as in the Hermitian case (that is with g = 0). 

2. If g < g < (jcr then a finite proportion of eigenvalues moves out of the real axes 
and places itself on very smooth curves in the complex plane. These curves converge to non- 
random limiting curves as n oo. Moreover we found the limiting density of states for non- 
real eigenvalues and proved that the "asymptotically" real eigenvalues have the same limiting 
distribution as in the case of the self-adjoint model. 

3. If g > gcr then virtually all eigenvalues leave the real axes. 

One thus concludes that the phenomenon described above persists as n is growing. The 
fact that the spectrum of the limiting operator is two-dimensional means that on this spectrum 
(but off the above mentioned curves) the resolvent of exists but its norm tends to infinity as 
n ^ oo. (See also [ p^ where the spectrum and the behavior of the norm of the resolvent of 
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random bidiagonal matrices were studied.) 

In the present paper we investigate spacings between the neighboring non-real eigenval- 
ues. We do this for models with selfaveraging (deterministicl) potentials which are introduced 
below. The class of selfaveraging potentials is very wide and includes in particular random 
stationary, quasi-periodic, and many other potentials. 

The approach used in [0, p was based on the theory of products of random matrices and 
on the potential theory. We would like to emphasize that the present paper is self-contained 
and that the technique we use here is much simpler and more straightforward than that we used 
before. It allows us: 

- To prove that the non-real eigenvalues of NHA model behave in a very regularly way: not 
only do they belong to very smooth curves as n cxo but also the difference between any two 
neighboring non-real eigenvalues z^+i and can be calculated with a great deal of precision as 

27ri 1 

Zk+l — Zk — —, r — , 

j[zk) n 

where f{z) is an analytic function of z (Theorem |33|). This is our principal new result. 

- To establish the log-Holder continuity of the density of states for Hermitian Anderson 
models with selfaveraging potentials. 

- To obtain (essentially as a by-product) much simpler proofs of our previous results listed 
above. 

Remarks. 1. In [|7|, ^ we considered tri-diagonal matrices with off-diagonal elements 
H(j,j + 1) and H{j,j — 1) depending on j. All results of this paper can be extended to 
these models. The main additional condition that is needed is the existence of finite limits 
lim„^oo YJj=i In H{j, j + 1) and lim„^oo n''^ YJj=2 ^{j, j - I). 

2. We have already mentioned above that, apart of the fact that stationary potentials provide 
natural examples of self-averaging potentials, randomness does not play any role in this paper. 
However as soon as one wishes to understand the real meaning of the next order term in formula 
( Plj ) randomness becomes crucial. The very same thing applies to the properties of asymptoti- 
cally real eigenvalues. Namely, it is natural to conjecture that the asymptotically real eigenvalue 
in fact are real for sufficiently large values of n. However there are reasons to believe that in 
order to prove this conjecture one should restrict himself to the class of random potentials with 
good properties. 

3. The approach based on the theory of products of random matrices (TPRM) is more 
difficult than that of the present work. However, all main results of this paper can be obtain 
within the framework of the TPRM. Moreover this is what we did first. The additional advantage 
of the TPRM approach is that the case of one-dimensional differential Schroedinger operators 
can be treated by this method in exactly the same way as the discrete case. The attempt to 
extend the approach of this paper to the case of continuous space would lead to the necessity of 
using a certain regularization procedure similar to the one used in [0]. 

The paper is organized as follows. We introduce the class of selfaveraging potentials and dis- 
cuss the log-Holder continuity of the relevant density of states in Section 2. Section 3 contains 
the statement of main results which are proved in Section 4. In Appendix A we prove a version 
of a well known property of Lyapunov exponents (which we use in Section 2). This allows us to 



3 



(a) make our exposition even more self-contained and (b) to demonstrate one more application 
of the technique we use. Appendix B contains an example of calculation of Lyapunov expo- 
nents for models whose potentials have rare high peaks. We believe that these examples are of 
their own interest but the initial purpose of finding them was to provide a natural completion 
for the study of selfaveraging potentials. 



2 The selfaveraging potentials 

We introduce now a class of deterministic potentials for which the distribution of the eigenvalues 
of the NHA model ([T|) - (0) will be investigated in this paper. 

Given an infinite sequence of real numbers q = {qk}'kLv consider the sequence of 
selfadjoint operators if °, n = 2,3, . . ., with potential q. Let Nn{E) be the distribution function 
of the eigenvalues of 

Nn{E) = : Ei G spectrum of and Ei < E}. (3) 

Definition We say that a real potential q is selfaveraging if the integrated density of states 
of the self-adjoint Anderson model with this potential exists, i.e. there exists a non-decreasing 
function N{E) such that Nn{E) — > N{E) as n ^ oo at the points of continuity of N{E). 

Remarks: 1. We have borrowed the name 'selfaveraging' from the theory of random op- 
erators, where it is normally associated with convergence of Nn{E) to a nonrandom limiting 
function. 

2. The class of selfaveraging potentials is very wide. For example, it contains decaying 
potentials, periodic and almost periodic potentials, and stationary random potentials^, see e.g. 
book [|T^ for proofs and more examples. 

3. We do not require J dN(E) = 1. However, in many cases (and in particular in those 
mentioned above) the sequence of measures dNn{E) is tight and hence cannot lose mass. 

4. The choice of periodic boundary conditions for is not essential, for example one 
could use the Dirichlet boundary conditions (po = (fn+i = instead. This is because changing 
boundary conditions amounts to a rank two perturbation of H^, and hence has no effect on 
N(E) provided the perturbed finite interval operators remain selfadjoint. Our preference for 
the periodic boundary conditions will become apparent in the next section. 

Let 

Unix,y)= / \n\x + ty-E\dNn{E) = -yin\x + ty-E,\, (4) 
J-oo nj^^ 

where the summation is over all eigenvalues of ii°, and 

/ + CXD 
ln\x + iy - E\dN{E), x,y eR, (5) 
'OO 

'if {qk}kLi is a strictly stationary sequence of random variables then Nn{E) is weakly converging for almost 
all realizations of q with respect to the corresponding probability measure. 
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Obviously the logarithmic integrals U{x, y) and ?7„(x, y) are the real parts of the analytic func- 
tions 

/ + 00 
\Yi{E-z)dN{E) = U{x,y)+iV{x,y). (6) 
-oo 

/ + 00 
\n{E - z)dN^{E) = Unix, y) + iVn{x, y). (7) 
-oo 

Here and below we consider analytic functions defined in the upper half -plane C+ = G C : 
Im^; > 0}, and the branch of the ln(£' — z) is chosen so that ln(— ?) = —if. 

The functions defined in (Q) - play an important role in this paper. In this section we 
study their properties under the following conditions 

CI The potential q is selfaveraging. 

C2 sup„>i i ELiln(l + \qk\)<C< +00. 

The main role of Condition C2 is to ensure that the functions U{x, y) and F[z) are well defined: 

Proposition 2.1 Assume C1-C2. Then for every x,y eM. the integral in is converging. 
Proof. Suppose first that y > and let z = x + iy. Note that 

Iny <Un{x,y) <ln{2+\z\) + C for all x G R and y > 0. (8) 
The LHS inequality is trivial, and the RHS inequality is ensured by Condition C2. Indeed, 

1 " 

Unix,y) < -^ln(2+ + \q,\), 
i=i 

where the last inequality follows, e.g., from the representation f/„(x, y) = ^\n\ det{zln — H^) \ 
and Hadamard's inequality for the determinants. By Condition C2, 



1 " 

-^ln(2+ |z| + < ln(2 + |z|) + C. 



n 
i=i 

If y4 > — OO and B < +oo are points of continuity of N{E) then, in view of Condition CI, 

nB pB i>B 

/ \n\z - E\dN{E) = lim / \n\z - E\dNn{E) < lim / \n\z + i - E\dNniE) 

J A "^^"^ J A J A 

< limsup 1/ + 1). 



Applying (|8|), we obtain that 

-B 



In \z - E\dN{E) <\n{3 + \z\) + C. (9) 



A 



As y 7^ 0, the function In |z — i?| is bounded from below, and (Q) implies that the integral in 
is converging for ?/ > 0. By the symmetry, it is also converging for y < 0. 

We shall now make use of the following inequality which will be proved later (Theorem 
U{x,y) > —Co for some cq > and all x and y ^ 0. This inequality together with the 
monotone convergence theorem yield that 

"+00 

(10) 



/+00 
\n\x - E\dN{E) > -cq, 
^00 



hence the integral in @ is also converging for y = 0. □ 

Remark. It follows from ( |T0[ ) that A^(-) is a continuous function. Therefore, under Con- 
ditions CI and C2, we have that Nn{E) converges pointwise to N(E) as n — > 00, and the 
convergence is uniform in E, —oo<E<oo. 

In view of the inequalities 

Qn-2<H'^<Qn + 2 (11) 

where Qn is the operator of multiplication by q, {Qn^)j = Qj^j, j = 1, ■ ■ ■ ,n. Condition C2 
also ensures that the sequence of measures dNn{E) is tight: 

Proposition 2.2 Suppose that q satisfies Condition C2. Then for any £ > there exists B > 
such that Nn{B) — Nn{—B) > 1 — efor all n. 

Proof. Denote by Xb{E) the indicator-function of the interval [B, +00), and let . . . , be 
the eigenvalues of if °. Then for any i? > 

1 . s 1 . ^ln(l + K■) 1^ , ln(l -f 2 -f g,) 

where the last inequality follows from ([TT|). We use here the following fact: if f{E) is a non- 
decreasing function and A < B then tr f{A) < tr f{B). Since X-b(2 + Qj) ln(l + 2 + qj) < 
ln(H- |2 + gj I) we conclude, in view of Condition C2, that Nn{B) > 1 — C'/ ln(l + i?) for some 
constant C' > and all n. Similarly, Nn{—B) < C" / ln(l + B) for some constant C" > and 
all n. □ 

We shall now investigate the relation between Fn{z) and F{z) in the limit n ^ 00. 

Proposition 2.3 Assume C1-C2. Then for every real x and y ^ 

liminf Un{x,y) > U{x,y). (12) 
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Proof. For all B large enough (so that [x — 1, x + 1] C [—B, B]) we have that for any n 

Un{x,y)> [ ln\x + iy - E\dNn{E) 

J-B 



and, by Condition CI, 



"B 

liminf f/„(x,?/) > / \n\x + iy - E\dN{E). 

-B 



By letting i? — > oo, we obtain ([12[). □ 

Under Conditions CI and C2, the sequence {f/„(x, y)} is not necessarily converging, even 
for y ^ 0, and for some selfaveraging potentials the inequality in ( |T2| ) is strict, for examples see 
Appendix 0. 

In view of (|8|), for any compact set K C C+, sup^ l-^n(-2)| < M{K) for some constant 
M{K) < +00 and all n. Hence the sequence {F„(z)} has a uniformly converging subsequence. 
We shall now describe all limit points of 

Theorem 2.4 Assume C1-C2. Suppose that Fn^{z) is a converging subsequence of {Fn{z)}. 
Then necessarily 

\iYa Fn^{z) = F{z) + c (13) 

for all z e C+ and some real constant c satisfying the inequality < c < In 3 + C, where C is 
the constant defined in Condition C2. 

Proof We note that = /+^(z-E)-MAr„(E) and F» = j^^{z~ E)-^dN{E). Since 

the function h{E) = (z — E)~^ decays to zero at infinity. Condition CI ensures that 

lim F:,{z) = F' {z) 

n— >oo 

uniformly in z on compact sets in C_|_. Therefore there exists a sequence of complex constants 
c„ such that Fn{z) — Cn — > F{z) as n — > oo for all z e C+. Passing on to the converging 
subsequence [z], we have that is also converging. Putting c = lim^ c„^ we arrive at (p^. 
It remains to prove that c is real, non-negative and satisfies the inequality c < In 3 + C. 

Due to our choice of the branch of the log-function we have that — tt < \m.\\i{E — z) < 
for all z E C+. Thus 

ImFniz) = J Imln(E - z)dNn{E) — > \mF{z), as n ^ oo, 

because the integrand is bounded, the sequence of measures dNn{E) is tight and we have Con- 
dition CI. Therefore the constant c is real. The fact that it is non-negative follows from Propo- 
sition P3| To complete the proof, note that Unj {x, y) converges to U (x, y) + cm the upper half 
of the xy-plane. Therefore, because of (||), ?7(x, y) + c < ln(2 + |z|) + C. Putting here x = 
and 1/ = 1, we obtain c < In 3 + c - f/(0, 1) < In 3 + C. (Obviously, U{x, 1) > for all x.) □ 

The following condition 
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C2* for any e > there is a i? > such that ^ ^ XbHijI) ln(l + \qj\) < e for all n 

guarantees the convergence of Fn{z) to -^(-z). It is obvious that Condition C2* is somewhat 
more restrictive than C2. On the other hand it is satisfied by many popular classes of potentials. 
For example, the random stationary potentials with finite expectation of ln(l + \qj) satisfy 
Condition C2* with probability one. It is also satisfied if 

1 " 

-^ln^+^(l + \qj\) <C* <+oo for some 5 > 0. 



Proposition 2.5 Assume CI and C2 *. Then 

Fn{z) — > F{z) uniformly in z on compact sets in C+ (14) 

n— +00 

and, in particular, 

Un{x,y) — > U{x,y) uniformly in z = x + iy on compact sets in (15) 



Proof. If the potential q is bounded then the statement of Proposition ^3] is a straightforward 
corollary of Condition CI and the fact that the functions Fn{z) are equicontinuous on any com- 
pact subset of C+. If q is unbounded then one needs to show additionally that the contribution 
of the tails of dNn{E) to Fn{z) is negligible in the limit n — > oo. Obviously, it will suffice to 
prove that: 

n 

for any £ > there is a 5 > such that ^J2XB{\Ek\) ln(l + \Ek\) < e for all n, (16) 

" k=l 



where the summation in dig ) is effectively over all eigenvalues of iJ„ such that \Ej\ > B. To 
complete the proof note that, in view of the inequalities in (PTj), it is apparent that Condition 
C2* implies (|T^. □ 

We finish this Section with a proof of the log-Holder continuity of N{E). This property is 
well known for random potentials and in this case it follows from the fact that U{x, y) > [0]. 
In turn, this inequality is a consequence of the Thouless formula according to which U (x, y) 
coincides with the Lyapunov exponent of with n = oo, see e.g. |T]]. 

In our case Conditions CI and C2 are too weak to guarantee the existence of the Lyapunov 
exponent even for non-real values of the spectral parameter, see examples in Appendix ^. How- 
ever these two conditions ensure that the function U (x, y) is bounded from below which, in turn, 
implies (very much in the same ways as in [Q]) that N{E) is log-Holder continuous. 

Theorem 2.6 Assume CI -C2. Then: 

( i) U {x, y) > —Co for some cq > and all real x and y. 
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(ii) N{E) is log-Holder continuous: for any E and \cr\ < ^ 

\N{E + a)-N{E)\=c{E,a)\ln\a\\~^ where lim c(E, a) = 0. (17) 



If E belongs to a compact set then c{E, a) < ci with the constant ci depending only on 
this compact set and the constant C in Condition C2. 

Proof Part (i). In view of (|T^ and the symmetry in y, it will suffice to prove the inequality for 
y > only. It follows from Theorem that 



liminf f/„(x,?/) = f/(a;,?/) + Co (18) 

n— ►oo 

for some cq > and all x and y ^ 0. To finish the proof, it is sufficient to show that the 
LHS in (|rS|) is non-negative. If the lim„^oo Un{x, y) exists then it coincides with the Lyapunov 
exponent, and hence is non-negative. The general case can be treated similarly (see Appendix 

Part (ii). Since / In |a; - E\dN{E) > -cq, we have that 

/ \n\x- E\ dN{E) + / hi\x - E\ dN{E) > -cq, 

J\x-E\<l J\x-E\>l 



and 

1 



In^ ^dN(E)< / \n\x-E\dN(E) + co<U(x,l)+co. 

\x - E\ J 

|x-£;|<l |3;-E|>1 

Therefore, for any \5\ < |, 

1 



Uix,l) + Co> / In^ -dN{E)>\N{x + 6)-N{x-6)\\n-^^, 

J \x-E\ \d\ 

\x-E\<\5\ 

and 

\N{x + S)-N{x-5)\ < [Uix,l) + Co]\\n\S\\'\ (19) 

Note that for any compact set C M, maxx U{x, 1) < +oo. This is because f/(x, 1) is 
continuous in x. 

Now, define for \6\ < ^ 

rx+S -1 

c(x,6)= / In^ —dN(E). 

Jx \x-E\ 

Obviously, 

|iV(x + 5)-iV(x)|<^. 

To complete the proof, note that ( |T9| ) implies that the measure dN{E) has no atoms, and there- 
fore 

c(x, 6) — > when 5^0. 

□ 
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3 Main results 



For the sake of convenience and clarity of exposition, we shall formulate and prove our results 
for the class of potentials q satisfying Conditions CI and C2*. We emphasize however that 
our main results hold true, modulo trivial modifications, under Conditions CI and C2, and the 
corresponding proofs are identical to those given in Section ^ This is a mere reflection of 
Theorem and the fact that our proofs are based on the convergence of Fn{z) to F{z). 



3.1 Notations and auxiliary statements 

Let us fix any finite interval [a, h] of the real axis {a < h). Most of our results apply to the part 
of the spectrum of iJ^ belonging to the strip {z : Re z G [a, h]} in the complex plane C. 
We define several critical values of the parameter g: 

g =inff/(x,0), (^cr = sup ?7(x, 0). (20) 

where we have introduced the notation S for the support of the measure dN{E), and 

9A0'M= inf U{x,Q), gcria,b) = sup U{x,0). (21) 

xesn[a,b] x&Sn[a,b] 



It may happen that gcr = +00, and it is obvious, in view of Theorem that _g^^ > for any 
potential satisfying Conditions CI and C2. 

For every g define 

Ag = {x eR: f/(x,0) < g}. (22) 

If fl' < fi' then Ag = 0, otherwise consists of (possibly infinitely many) disjoint open 
intervals: 

j 

We note also that U {x, y) = U (x, —y) and that 

d 

— U(x,y)>0 for any a; G M, y > 0. (24) 

dy 



Proposition 3.1 Suppose that g > g^^ and let {aj, bj) be the intervals defined in (|23|j. Then the 
level set Cg = {{x,y) : U (x, y) = g, y > 0} consists of disjoint analytic arcs 

y = Vjix), aj <x < bj, (25) 

whose end-points lie on the real axis, i.e. yj{aj + 0) = yj{bj — 0) = 0, if —00 < aj, bj < +00. 
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Proof. If xo ^ Ag then U{xo,y) > g for all 7^ 0. Therefore the equation U{xo,y) = g cannot 
be solved for y > 0. 

Consider now any of the intervals {aj, bj). If xq E bj) then U (xq, 0) < g, and in view 

of and U{x,+oo) = +00, there exists a unique positive solution yo *= Vji^o) > of 
the equation U{xo,y) = g. As U{x,y) can be analytically continued into a neighborhood of 
(xo, 1/0) in C^, the implicit function theorem asserts that yj (x) is analytic in a disk \x—xo | < 5 in 
the complex x-plane. The union of all such disks, when xq runs through {aj, bj) covers (aj, bj). 
Therefore the function yj (x) can be analytically continued into a domain in the complex x-plane 
that contains (aj, bj), and, for any closed interval [a, /3] C (aj, bj), this domain contains 

Da^p = {xeC: a - h <Rex < (3 + h, \lmx\ < h} (26) 

for some h > 0. 

If aj > —00 then yj{aj + 0) = 0. For, if not then y := Um sup^.^^^ ^.^^^ Vji^) > 0. But 
then U{aj,y) = g and hence U{x,0) < g for every x from some neighborhood of aj which 
contradicts the definition of aj as the end point of our interval. The same argument proves that 
if bj < +00 then yj{bj - 0) = 0. □ 

3.2 Statement of results 

We are now in a position to formulate our main results. 

Theorem 3.2 For any g > all the eigenvalues of belong to the level lines of the function 
Un{x, y) defined by the equation 

Un{x,y)=g + -\Yi{l-e-^^). (27) 



Theorem 3.3 (i) Suppose that g < g^^ici, b). Then for any £ > there exists uq = no{e, g, q, a, b) 
such that for any n > uq all the eigenvalues Zj of with Re Zj G [a, b] belong to the e- 
neighborhood of the real axis: \ l^azjl < e. 

(ii) Suppose that g > and [aj, bj) is one of the intervals comprising Ag. Then for any 
[a, P] C {aj, bj) there exists rii = ni{q, g, a, (3) such that for any n > rii there exists a solution 
yj^n{x) to equation ( ]2^ which is analytic in the domain Da^p defined in ( |2^ ) and 

lim yj,n{x) = yj{x), uniformly in x G Da^p. (28) 

The function yj^n{x), for n > ni, is the only solution of^^ which is non-negative when x G 
[a, (3]. 

Remarks. 1 . The previous two theorems imply that if has eigenvalues in the strip aj < 
a < Re z < P < bj, then, for n> rii, they must lie on the analytic arc 

An{a, [3) = {{x,y) : y = %-,„(x), a < x < (3} 
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and on its reflection with respect to the real axis. 

2. Relation (Q) implies that the arcs An{cy, (3) converge to the level lines of U{x, y) when 
n ^ oo together with all their derivatives. 

3. We did not make use of two out of the four critical values introduced in ( ^0[ ) and (pT]). 



However, their role is clear: if g > 5'cr(a, b) then no limiting eigenvalue curves grow out of the 
support of 5 n [a, b]. If (ja is finite then [a, b] can be replaced in this statement by S. 

The next two theorems describe the asymptotic distribution of the eigenvalues of along 
the arc An{(y,P). In particular they state that H^, for large n, does have eigenvalues on 

By P) we denote the number of complex eigenvalues of H^^ lying on Anicn, (3). 
Theorem 3.4 For any closed interval [a, [3] C (a^, bj), 

lim -u^{a,(3) = ^[e{(3)-e{a)l 
n^oo n 2n 

where 9{x) = —V{x, yj{x)) and V{x, y) is the imaginary part of F{z). 

Remark. Let / be the natural parameter on the curve y = yj{x), that is the length of the part 
of this curve contained between say (a,yj{a)) and {x,yj{x)). A simple calculation involving 
the Cauchy-Riemann equations for F{z) shows that d6 = \f{z{l)\dl, where 



Hence 



1 1 /"^ 1 /-z^+^wC/^) 

lim -Ma,P) = — / 9'{x)dx = — / (30) 



la+iyj{a) 

where the integration is carried out along the path y = yj{x) from a + iyj{a) to P + iyj{P) 



Theorems and are not entirely new and can be inferred from Theorems 2.1 and 2.2 
in [^. We are now going to formulate our principal new result. Let [a,P] be the same as 
before. Let us label the eigenvalues = Xk + iyt of lying on the arc so that 

OL < Xi < X2 ■ ■ ■ < Xm P (we note that in fact the multiplicity of these eigenvalues is one 
and the inequalities here are strict; this follows from the inequality 6^{x) > Cq > which is a 



part of the proof of Theorem |J!4| ) 



Theorem 3.5 For any two consecutive eigenvalues Zk and Zk+i of lying on An{oi, P), 

n{zk+i - Zk) = + Sn{zk, Zk+i) (31) 

where 

lim 6n{zk,Zk+i) = uniformly in Zk,Zk+i G An{a,P). (32) 
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4 Proofs 



The eigenvalues and eigenfunctions of are determined by the equation 

-e'^ipk+i + QkVk - e'^ifk-i = zipk, l<k<n (33) 

were 

V^o = V^n, (fil = (fin+l- (34) 



The parameter g can be eliminated from ( P3| ) by making use of the standard substitution ifk 
e'^^-ipk which transforms (Pl) into 



-^pk+i + Qk^k - ^k-i = z^pk, I <k <n, (35) 
and boundary conditions ( |3^ into 

V^o = e-^^V-n, ^1 = e-"^7/>„+i. (36) 

Note that the transformed boundary conditions are asymmetric (unless g = 0). 

To solve equation ( 05] ) we shall follow the standard routine and rewrite it in the matrix form: 



Then 



On the other hand, 



^n+l\ _ no 



because of boundary conditions (^6|). Therefore the eigenvalues of are determined by the 
equation 

det\Sn{z)-e''n] = Q. (37) 
Since det Sn{z) = 1 for all z, we have that 

det[Sn{z) - e"3/] = 1 - e'^Hr Sn{z) + e^""^. 

Hence: 

Lemma 4.1 2; is an eigenvalue of iff 

tr5„(z) = e"^ + e-"^. 

The trace of the matrix Sn{z) is a polynomial in z of degree n. The following representation 
of this polynomial, which is well known in the context of the discrete Hill equation (see e.g. 



[ ]13| ] or [^) is useful for our purposes. 
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Lemma 4.2 Let Ej, j = 1, . . . ,n, be the eigenvalues ofH^. Then 



iiSn{z) = \{{Ej- z)+2. (38) 



Proof. For (7 = 0, Lemma pTTI asserts that the polynomials nj=i(-^i ~ ^) ^'^^ 'S'„(2;) — 2 have 
the same set of zeros. It is easy to verify both polynomials have the same coefficient, (—1)", in 
front of the highest power of z, hence they must coincide. □ 

Here is our main technical lemma: 



Lemma 4.3 Suppose that g > 0. Then z is an eigenvalue ofH^ iff 



, , 2 / Ii£ -"9 \ tTT 

FJz) = — In 62—6 2 H 

n \ / n 



mod 



2Txi 



n 



where Fn{z) is the function defined in (0). 

Proof. It follows from Lemmas ^]T] and that z is an eigenvalue of iff 



(39) 



- z] 



ng -ng r, 

(62-62). 



Since Fn{z) = ^ Yl%i — z), equation ( |40| ) is equivalent to (^, provided g 0. 

We are now in a position to prove Theorems and |T3| . 

Proof of Theorem ^Jj This theorem is a straightforward corollary of Lemma p3t 

Proof of Theorem P3| : (i) Let (7(e) = min^gja f,] U{x,e). According to Proposition |23 
find 77-0 such that for all x E [a, 6] 



(40) 



□ 



one can 



\Unix,e) ~U{x,e)\ < -[g{e) - g^^{a,b)] 



if n > (Note that g{e) > g^^ia, b) because of (|24l).) Thus, for all a < x < 6 and y > e, 

Un{x,y) > Unix,e) > U{x,e) - ^[^(e) - _g^^{a,b)] > ^U{x,e) + ^g_^^{a,b) > _g^^{a,b). 

Recall that by the assumption g^^ia, b) > g. Since g > g + ^\n{l — 6""^) for any n > 0, we 
conclude that, for all n > uq equation ( p7| ) has no solutions in the half-strip a < x < b, y > e. 
To complete the proof remember that f/„(x, —y) = Un{x, y). 

(ii) First, note that for every real x equation ( p7| ) has one non-negative solution at most. 

Now, let g > g^^ and [a, P] C bj) where (aj, bj) is one of the intervals comprising A^. 
Recall that yj{x) is analytic in Da,f3, see Proposition |3.1|. Because of the compactness of [a, P], 
it will suffice to prove the existence of the solution yj,n{x) to equation (^71 ) and its convergence 
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to yj{x) as n — > oo in a small neighborhood of every point {x,yj{x)) where x runs through 

Fix X G [a, P] and consider the point (x, y) where y = yj{x). It follows from the integral 
representations for Un{x, y) and U (x, y) that these two functions are analytic in the domain 

D^={{x,y): < |, < |}. 

We shall use the following general lemma. Put =^ {(x, y) : \x — x\ < r, \y — y\ < r}. 



Lemma 4.4 Let ^{x,y) and ^{x,y) be two functions analytic in D,. and such that for all 

{x,y)eDr 

\Kix^y)\ < ci, < C2 < \^'y{x,y)\ < C3, \^{x,y)\ < 1. 

Suppose that y) = 0. Then there is a positive Eq which depends only on ci, C2, C3, and r 
(but not on $(-, ■), $(■, •)/) such that the equation 



^{x,y) +e^{x,y) = 



(41) 



has a unique solution y = y{x, e) which is analytic in (x, e) in the domain {(x, e) : |x — x| < 
2eo, \e\ < 2eo} and y{x, 0) = y. 



Proof of Lemma Consider the function G{x, e, y) =^ y) + ^^(x, y) of three complex 
variables x, e, and y. In the domain Dz we have: 



<i:.(x,y) =(27r)-i 



\u\=r 



< 



and similarly |$^| < ^. Hence, for e sufficiently small, G'.j. and G'y are close to and $^ 
correspondingly. It is clear that | | < 1 . The implicit function theorem for an analytic function 
of three variables x, e, and y implies now the existence of the solution y = y{x,e) to the 
equation G{x, e, y) = 0. It should be emphasized that the domain where this solution exists and 
is analytic depends only on the corresponding estimates of G'^,G'y, and G'^. □ 

To finish our proof of Theorem |^, note that in our case U (x, y) plays the role of <l>(x, y) 
and equation (|4l|) has the form 



U{x,y)+e 



Un{x,y) -U{x,y) 



0. 



^0 



Here we first choose Eq so that to satisfy the conditions of Lemma and then choose uq such 
that \Un{x, y) — f/(x, y)\ < Sq for all (x, y) e D and n > uq. The wanted result follows from 
our Lemma when e = Eq. □ 



Define 



On{x) = -Vn{x,yj^n{x)) and 9{x) = -V{x,yj{x)) 
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for X G [a, (3\ C (flj, hj) and n > ni with ni as in part (ii) of Theorem As before, Vn{x, y) 
and y) are the imaginary parts of the analytic functions Fn{z and F{z), see (@) and In 
view of Theorem and Proposition |XT], we have that 



hm 6n{x) = 6{x) uniformly in X G [a, [3]. 

?i— >oo 

It follows from the Cauchy-Riemann equations for and F{z) that 



(42) 



|Vf/„(x,y)p 



9j/ 



and 6 (x) 



\VU{x,y)\' 



^U{x,y) 



dy 



(43) 



Therefore we also have that 

lim 9'^{x) = 9'{x) uniformly in a; G [a, P]. 



(44) 



As -^Un{x, y) and ^U{x, y) are positive in the upper half of the xy-plane, the functions 9n{x) 
and 9(x) are monotone increasing. Moreover, in view of (P^), it apparent that there is a constant 
Co > such that 

9'n{x) > Co > for every x G [a, (3] and n > ni. (45) 
Now we are in a position to prove Theorems ^j^and ^31 . 

Proof of Theorem On the arc An{a,(3), i.e. when z = x + iyj^n{x), a < x < (3, the 
eigenvalue equation ( l39|) reduces to 



Jndn (x) 



-1. 



(46) 



When X runs through [a,P] in the positive direction, 9n{x) gets positive increment and w 



moves anticlockwise along the unit circle |w| = 1. Obviously, (3), the number of 
eigenvalues of on (3), is, up to ±1, equal to the number of circuits completed by w 

when X completes its run. Thus 



Vn{0l,(3) 



n[9n{oi) - 9n 
27r 



where \k\ < 1. When n oo, 9n{x) converges to 9{x) = —V{x,yj(x)), and, therefore, 
iz/„(a,/3) converges to ^[9{(3) - 9{a)]. □ 



Proof of Theorem p3t Let zi = xi + iyj{xi) and zi^i = xi+i + iyj{xi+i) be two consecutive 
eigenvalues of on Ania, (3). We assume that x^+i > x^. It follows from equation (^^ that 



6'„(a;i+i) - 9n{xi) 



n 



and therefore 



27r 



n 9'^{x*) 



(47) 
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for some x* G (x/, xi+i). In view of (P3|), 



27r 1 

< -xi< — — . 

Con 

for all n large. Hence x* ^ x/ as n — > oo, and ( P7| ) and (Q) imply that 

27r 



n(xi+i - xi) 



e'{xi) 



+ 5„(xz,Xi+i) 



where 



lim (5„(xi,x;+i) = uniformly in Xi, Xi+i G [«,/?]. 
To prove - @, note that 

Zl+l - Zi = X;+i - X; + Z?/^- „(x**) (x^+i - X;) 

for some x** G (x;, x^+i). By making use of (|^), one obtains that 

n(2;i+i - Zi) = ^^[1 + Wj,n{^**)] + ^n{xi, Xl+i). 

Now ( prt ) - easily follow from Theorem ^3] and the following relation 



l + iyJx) 



1 



e'{x) iF'{z) 
To verify this relation, make use of the equation 



z=x+iyj{x) 



j \og\x + iy,{x) - E\dN{E) = g, 



to obtain 



Now, in view of (|43|), 



1 +iyAx) 



1 



e'{x) Uy{x,y) +iU'^{x,y) 
and ( |50| ) follows by the Cauchy-Riemann equations for 
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A Appendix 



Proposition A.l For all real x and y ^ we have 

liminf f/„(x, y) > 



(51) 



Proof. This result can be proved in many ways. We present here a proof based on (p8|). 
It follows from (pb that 



1 



1 1 

-lntr5'„(z) - Fn{z) = - hi 
n n 

Therefore for every z E C+ such that Imz > 1 

1 



n-=i(^.-^) 



(52) 



hm 

n— >oo 



n 



IntiSniz) - Fn(z) 



0. 



(53) 



The two functions in the LHS in ( p2| ) are analytic and uniformly bounded in n on compact sets 
in C+. Therefore, by the Vitali theorem, (Q) must hold for every z E C+. 

Consider now the eigenvalue equation for Sn{z). If z then Sniz) has no eigenvalues on 
the unit circle. As det Sn{z) = 1, we then have that for every non-real z the 2x2 matrix Sniz) 
has one eigenvalue, Xniz), in the exterior of the unit circle, i.e. \\n{z) \ > 1, and the other one, 
l/\n{z), in the interior of the unit circle. Thus tr Sn{z) = Xniz) + \n^{z) and 



— liaXnU) In tr Sn (z) 

n n 



n 



ln[l + A-^(^)]. 



(54) 



It follows from ( p8| ) that | trS'„(2;)| grows exponentially fast with n provided | Imz| > 1, and 
then so does the dominant eigenvalue of Sn{z). This is because | tr ^^(^r)! < |A„(2;)| + 1 (In 
fact, \n{z) grows exponentially fast with n for every non-real z.) Hence, for every z E C+ such 
that Imz > 1, 



lim 

n— >oo 



In XJz) 



n 



n 



0. 



(55) 



The two functions in the LHS in ( |54[ ) are analytic and uniformly bounded in n on compact sets 
in C+. Therefore, by the Vitali theorem, ( |55| ) holds for every z E C-j.. From ( |53| ) and ( |55| ) we 
have that 

-\n\n{z) - Fn{z) =0 



lim 

n— »oo 



n 



for every z E C+. Taking the real part 

1 



lim 

n— >oo 



n 



In |A„(a; + iy) \ - Un{x,y) 



0, 



and therefore (recall that |A„(2;)| > 1) 



liminf Un{x, y) = liminf — In |A„(x + iy)\ > 0. 

n— >oo n— >oo n 



□ 
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B Appendix 

Obviously the integrated density of states, N{E), depends on the potential q. To make this 
dependence explicit, we shall write in this section Nn{E; q) and N{E] q) instead of Nn{E) and 
N{E). 

Let fci, /c2, • • • be an increasing (infinite) sequence of natural numbers such that 

#{j : k, < n} 

> as n ^ oo, (56) 

n 

and let v = {f^}^^ be a potential supported by the sequence kj, i.e. = unless k E {kj}. 

Proposition B.l Ifq is a selfaveraging potential then so is q = q + v, and N{E] q) = N{E, q). 

Proof. According to the well known theorem from linear algebra, if A and B are two self adjoint 
n X n matrices then the number of eigenvalues of the matrix A + B in interval A differs from 
that of the matrix A by rank(A — B) at most. Hence 



\Nn{E-q)-N4E-q)\ < 



#{j : k, < n} 



n 

which proves the proposition. □ 

Let us now assume that the potential v satisfies Condition C2 and |f oo as j oo. 
Define ^ 

Sn = - ln(l + \VkA). 



n 

j:kj<n 



and 

r" + 00 



Un{x,y)= ln\x + iy - E\dN{E]q), y > 0. 



Theorem B.2 Let q be a selfaveraging potential with sup^ \qk\ = M, and q = q + v. Then 

lim [Unix,y) - = U{x,y). 

n— >oo 

Proof. We prove this Theorem under the following additional condition on v: |f fcj_|_i | > l^fc^ | + 
2 + M for all j. The proof for the general case requires minor but cumbersome modifications. 
Let z = X + iy.By definition, 

1 " - 1 ~ 1 

Un{x,y) = -y^\n\z - Ek\ = - \n\z-Ek\ + - \ii\z - Ek\, 

^==1 |^fc|<2+Af \Ek\>2+M 

where the E^ are the eigenvalues of with the potential q. By Proposition JB.l| the first sum 
converges to U{x,y) when n ^ oo. We note next that the eigenvalues Ek in the interval 
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\Ek\ > 2 + M have the following property: for all but may be a finite number of them there 
is a unique j such that \Ei: — v^ ] < 2 + M. This is due to the fact that the eigenvalues of the 
operator of multiplication by v differ from the eigenvalues of with the potential g by 2 + M 



at most. Hence taking into account that In j 



'-n^ when \Ek 



oo, we obtain that 



lim 

n^oo 



\E,,\>2+M 



kj<n, \vk- \>2+M 



It is apparent that 



lim 

n— >oo 



- \n{l + \VkA) - Sn 

n ^-^ 



kj<n, , \>2+M 



which completes the proof. 



□ 



It is easy now to construct examples showing that the statement of Proposition cannot 
be improved. 



Example 1. Let kj = and Vj2 = , j = 1,2,.... Then lim 



lim^^oo Unix, y) = U{x,y) + \. 



Example 2. Let kj = 2^ and = j = 1,2, Then limsup^^g^s 



n— >oo '^n 



T. Hence 



2 and 



liminf„^ooS„ = 1. Hence lim sup„^^ t/„(a;, = f/(a;, + 2 and liminf^^oo Z/) = 
U{x,y) + 1. It is easy to check that for every 1 < c < 2 there is a subsequence Unj{x,y) c! 
converging to U (x, y) + c when j ^ oo. 
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